Charge Pumping and Photovoltaic effect in Open Quantum Dots. 
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We propose a random matrix theory to describe the influence of a time-dependent external field 
on electron transport through open quantum dots. We describe the generation of the current by an 
oscillating field for the dot, connected to two leads with equal chemical potentials. For low frequency 
fields our results correspond to adiabatic charge pumping. Finite current can be produced if the 
system goes along a closed loop in parameter space, which covers a finite area. At high frequency a 
finite current is produced even if the loop is a line in parameter space. This result can be explained 
in the same way as adiabatic pumping but considering the evolution of the system in phase space 
rather than in parametric space. 

PACS numbers: 73.23.Ad, 72.15.Rn, 72.70.+m 



I. INTRODUCTION 

Adiabatic charge pumping through open quantum dots 
was studied recently in the literature both theoretically 
and experimentally j|. Such pumping occurs in a 
system described by a Hamiltonian periodic in time with 
a period T p larger than all other characteristic time scales 
of the system. After one period, the system returns to 
its initial form; however charge Q can be transmitted 
through a cross-section of the system: 



Q — IdcT p — 



(I(t))dt, 



(1) 



where (...) denotes quantum mechanical and thermody- 
namic averaging. 

To obtain a finite transmitted charge at low frequen- 
cies, the Hamiltonian should depend on at least two pa- 
rameters. In Refs. fl]-||] the time dependence of the 
Hamiltonian was replaced by a dependence on param- 
eters and the system was considered quasistationary for 
each parameter value. The transported charge during one 
period of the Hamiltonian was calculated as an integral 
in the parameter space. The theory || shed some light 
on the recent experiments 0, namely on the amplitude 
dependence of the root mean square fluctuations of the 
transmitted charge, averaged over different realizations 
of the Hamiltonian. 

A very similar phenomenon was considered previously 
by Falko and Khmelnitskii, who theoretically studied the 
photovoltaic effect in mesoscopic micro junctions |q]. The 
experimental observation is described in Ref. 0]. The 
photovoltaic effect is a generation of d.c. -current by radi- 
ation of a finite frequency. (It is obvious that this effect 
can only be non linear in the oscillating field.) The bi- 
linear regime of adiabatic pumping, is precisely the 
circular photovoltaic effect introduced in Ref. M and ap- 
plied to a mesoscopic system in Ref. [pi. The results of 



Ref. H are not directly applicable to quantum dots be- 
cause in micro junctions the Thouless energy Et ~ 1/T erg 
is of the same order as the inverse escape time l/r eS c, 
whereas for quantum dots 1 / r csc <C -Et • (Here r org is the 
characteristic time for a classical particle to cover all of 
the available phase space in the dot and we put h = 1.) 
Therefore the considerations of Refs. [QJ|] have their own 
physical significance. On the other hand, the theory of 
Ref. || is not restricted to the adiabatic regime, the re- 
sults being valued in a broad interval of frequencies. 

The purpose of the present work is to go beyond the 
adiabatic approximation for d.c. -current generation in 
open quantum dots. One can identify two contributions 
to the d.c. -current — reversible and irreversible. To make 
a connection with the terminology of the photovoltaic 
effect, used in Ref. Hj, we consider the bilinear d.c- 
current response through the dot, generated by several 
time-dependent perturbations ifi(t) = ip% u e . 

The direct current Tdc can be written at to — > as 

Joe = iueijipi^V*j,u + u 2 jij(p itU} cp*j jU , + C(w 3 ), (2) 

where ey and jij are real sample specific coefficients. 
They are not fixed by any symmetry (a sample does not 
have any!) except the condition that 7dc is a real num- 
ber, which gives the requirements 



(3) 



These relations mean that as the direction of the con- 
tour in the parameter space {(fi(t)} is reversed (<£i(t) — > 
tpi(-t) or ipi tbJ — > <p*u)i t ne nrs t term changes its sign 
whereas the second term remains intact. In the language 
of the photovoltaic effect, the first term is the circular 
photovoltaic effect, and the second term is the linear ef- 
fect, fri Equation makes an explicit connection be- 
tween the bilinear responses of d.c.-current generation 
through open quantum dots and the photovoltaic effect. 
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We will see that the separation onto reversible and 
irreversible parts goes far beyond the bilinear and low 
frequency expansions. We will find that the reversible 
part vanishes at high frequency whereas the irreversible 
part saturates. We will also show that the irreversible 
contribution may be interpreted in a manner similar to 
adiabatic pumping. In adiabatic pumping the transmit- 
ted charge was determined by a contour in parameter 
space. The irreversible current is determined by the con- 
tour in the extended phase space {y>i{t),<fi(t), . . .}, see 
Sec. V. 

The first term in Eq. (^) vanishes for a single pump 
due to the antisymmetry of e^-, Eq. (||). The contour in 
parameter space degenerates to a line in this case, while 
in phase space the contour encompasses a finite area. 
We will see that the current is proportional to this area. 
Note, that the contour is invariant with respect to time 
inversion. 

The remainder of the paper is organized as follows. In 
section II we formulate the model which is studied in 
the paper. We mainly discuss the high temperature limit 
T 3> 10, where uj is the frequency of the perturbation. 
In section III we calculate the ensemble averaged fluc- 
tuations of the current. As an example, low frequency 
asymptotic of reversible (adiabatic) and irreversible cur- 
rents are considered. In Sec. IV we discuss how the 
irreversible contribution can be represented in the form 
of a linear integral over the contour in extended phase 
space. In Sec. V we discuss the low temperature limit 
and show how to take heating into account. Section VI 
summarizes our results. 



II. THE MODEL 

We consider the following experimental realization of 
the model. Gates near a two-dimensional electron gas 
(2DEG) form the shape of the dot. An oscillating volt- 
age is applied to the gates V\ (t) and V2 (t) . As a result of 
motion of electron energy levels in the dot, a current flows 
through the dot. The direction of the current depends 
on the particular realization of the dot and is zero on 
average. We calculate the fluctuations of the d.c. -current 
with respect to different realizations of the dot. 

Calculations will be performed for an open quantum 
dot in the limit of a large number of open channels A c h 
connecting the dot to the leads. This condition allows us 
to neglect the electron-electron interaction, which gives 
corrections of the l/N^ h order, (See Ref. ||]). The same 
condition permits the use of a diagrammatic technique, 
similar to that described in || , to calculate ensemble av- 
eraging. 

We also assume that the quantum dot is small and the 
Thouless energy Et ~ l/ r org is much greater than all 
other energy scales of the problem. In this limit, one 



can use random matrix theory (RMT) to study trans- 
port and thermodynamic properties of the system, see 
Ref. pp| . All corrections to the RMT are as small as 
Ach/ffdot, where (/dot = Et/S\, S± is the mean level spac- 
ing. We neglect the fluctuations of the time dependent 
perturbation from sample to sample, created by the volt- 
ages Vi^(t), since they depend on the small parameter 

<?dot « 1 - 




T v 2 (o 

FIG. 1. An experimental setup. The voltage applied to 
the gates Vi(t) and Vi{t) changes the shape of the dot, re- 
sulting in motion of the energy levels of the electrons in the 
dot. 

The Hamiltonian of the system is |lC| 

H = H D + H L + H LD , (4) 

where H^, is the Hamiltonian of the electrons in the dot, 
which is determined by the M x M matrix H nm : 

M 

Hb = ^ ^nHnmlpm, (5) 
ro,m=l 

where the thermodynamic limit M — > 00 is assumed. 
The coupling between the dot and the leads is 

#LD = (Wn^UWn + H.C.) , (6) 

where ip n corresponds to en electron state in the dot, 
tpa(k) denotes electron state (a, k) in the leads, and k 
labels the continuum of momentum states in each chan- 
nel a. For a dot connected by two leads with N\ and N T 
channels respectively, we denote the left lead channels by 
1 < a < N\ and the right channels by N\ + 1 < a < A c h, 
where A c h = N\ + N T . The coupling constants, W na , in 
Eq.© are @: 

rrr _f^n, H U = 01 < N<&, , s 

™ a ~\0, otherwise [ '> 

and r„ are defined below in Eq. ( ]f9| ) 

The electron spectrum in the leads near Fermi surface 
can be linearized 

H L = v F J2Wl{k)Mk), (8) 

a . k 
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where vp = 1/2-kv is the Fermi velocity, and v is the 
density of states per channel at the Fermi surface. 

The current through the dot is given in terms of the 
scattering matrices S(t, t') by the following expression, 
see appendix A: 



(/(*)) =e^A QQ Jdhdt-2 



(9) 



s «e (*. *i)//s (*i - *a)3L (*a , *) - /a (+»0) 

where / Q (i) is the Fourier transform of the electron dis- 
tribution function in the leads, (. . .) stands for the quan- 
tum mechanical and thermodynamic averages for a given 
ensemble realization (no ensemble averaging!) and 



A a /3 = S a f 




if 1 < a < Nr, 
if Ni < a < N ch . 



(10) 



The scattering matrix of the system, S(t, t'), is 
S af3 (t,t') = 5 aP 5(t - t') - 2mvW^GW(tJ)W m p, (11) 
and the Green's function G^l(t, t') is the solution to: 

i*M ~ ^ (<) + l7r "WW^ & R) (t, t') = S(t - t'), (12) 

where the matrices H and W are comprised by their ele- 
ments (see Eqs. (||) and @))- An analogue of Eq.(0) was 
used before in the energy representation in Ref. pH . 

Below we consider the special case of electrons in both 
leads being described by identical distribution functions 
fj(t) = f(t). The function f(t) is the Fourier transform 
of the Fermi-Dirac distribution function: 



+oo 



>T 



2 sinh nTt 



(13) 



In this case, we can derive (see Appendix B) another 
formula for the current through the dot: 



(I(t)) = 2eiirv r Tr / Jdtidt 2 f{ti - t 2 ) (14) 

x {w^GW(t,ti) Lff(ti) - H(t 2 )] G (A \t 2 ,t)WAj, 

which is more convenient for further calculations. 

We calculate the variance (mean square) of the trans- 
ported charge Q through the dot. We assume, that the 
Hamiltonian of the dot H in Eq.(|5|) is represented by a 
time-dependent matrix in the form: 



H(t)=H + V 1 cp(t) + V 2 ip(t). 



(15) 



Here, the time independent part of the Hamiltonian i~L is 
a random MxM matrix, which obeys Gaussian statistics 
with the correlator 



(16) 



where (...) means averaging over ensemble of matrices 
H, A = M{5 1 /tt) 2 and A' = A(l - g h /4M), and g h de- 
fines the crossover from orthogonal (g^ = 0) to unitary 
(.9h = 4M) ensembles. The parameter g^ has the mean- 
ing of the dephasing rate due to an external magnetic 
field in units of the level spacing Si, ]Io|Jl^| . It can be 
estimated as <?h ~ g^ot (^/$o) where $ is the magnetic 
flux through the dot and $ = hc/e is the flux quan- 
tum. The time dependent perturbation is described by 
symmetric M x M matrices, Vi, and the time dependent 
functions tp(t) and ijj(t). 

We make the simplifying assumption that the matri- 
ces V are traceless TrV — 0. In this case, the unknown 
parameters CV, are defined as 



M 2 8 x 



(17) 



where we used the fact that the matrix V is symmetric. 
(We note that according to the definition of Eq.(|l7|), C 
has dimension of energy.) The parameters C are also 
related to the typical value of the level velocities which 
characterize the evolution of an energy level e v (X) under 



the external perturbation XV, 1 1 3 1 



2ji _ de v dt v 

~V 13 ~ dXi dx 



, dXi dXj 



(18) 



Since all other responses (e.g., parametric dependence 
of the conductance of the dot) are expressed in terms of 
universal functions of the same parameters Cy [fl3| , they 
can be found from independent measurements. Calcula- 
tions can be generalized to the case in which the trace of 
the perturbation matrix is not zero. See, e.g. ||. 

For reflectionless contacts, the ensemble averaged scat- 
tering matrix S a p is zero. This condition determines the 
factors T„ in Eq. (Q): 



r„ = 



MSi 



(19) 



The ensemble average of the transmitted charge Q, de- 
fined by Eq. ([!]), is zero: Q = 0. We calculate the second 
order correlator with respect to an ensemble of random 
matrices, Eq. (pf): 



T D T D 



Q 2 



dtdt'{I{t)){I{t')). 



(20) 



For this purpose we use a diagrammatic technique, which 
has been applied to a similar class of problems. See Refs. 
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HQ. Here we present the basic elements of the dia- 
grams which will appear in sections III and V. The en- 
semble averaged Green's function G^m m the absence of 
perturbations is equal to 



G n m (e) — i 




iV ch ± it 
1 + , 7 , N ch <n<M: 



AM 



l<n< N ch . 



(21) 



Above we introduced the dimcnsionless energy, e, mea- 
sured in units of ^/A/4M = Si/2tt. We expanded these 
Green's functions in e/M and N^/M, since only those 
terms survive the thermodynamic limit M — > oo. For 
the same reason, in the expression for Gn with n < N c h 
one has to neglect such terms since the contribution of 
these elements to the final result is already of the order 
of N ch /M. 

The other element of the diagram technique used in 
the paper is an amputated average of the products of 
two Green's functions (see Fig. 2(b)), called the diffuson 
T>(ty, t2, t) and defined by 



(Gnm(t^ ,t 2 )Gmn(t 2 ,ty )) amp 

ty + ty f 2 + t 2 



(22) 



4MXV 



x 5{t+-t+-t-+t-) 



v - G 

, Uy i 2 



We can use this relation since the time arguments of diffu- 
son satisfy ty — t 2 = ty — t 2 . Introducing new variables 
<i.2 = {ty 2 + ty 2 )/^ an d t = ty 2 — ty 2 we obtain the 
following equation for the diffuson V(ty,t2,r): 



Wy +jCv{h ' T 



V(ty,h,T) = 6{ty-t 2 ), 

K v (t, t) = N ch + $ T (i, T)C${t, r), 



where 

*(i,r) = 

The solution to the above equation is 



*+§W(*-5) 



(23) 



(24) 



V(ty,t 2l T) = 9(ti -t 2 )exp 



K.v{t,T)dt . (25) 



Equations (|2^)-(p5|) are written in dimensionless vari- 
ables, so that energy and time are measured in units of 
Sy/2TT and 2it/5y respectively. Below, intermediate ex- 
pressions are also written in terms of dimensionless en- 
ergy and time variables, while the final answers are rep- 
resented in terms of physical quantities. 
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FIG. 2. (a) Diagrams for the ensemble averaged Green's 
function. The second term in the self-energy includes an in- 
tersection of dashed lines and is as small as l/M. (b) The 
Dyson type equation for the diffuson, "D(t~t , t j~, t% , t%). 




t 2 

— =f(t l -t 2 )[H(t l )-H(t 2 )] 
ti 



FIG. 3. The diagram representing the contribution to the 
charge correlator, Q 2 at high temperature T. 



III. MESOSCOPIC FLUCTUATIONS OF THE 
CURRENT AT HIGH TEMPERATURES 

In this section we will consider the high temperature 
limit, in which the frequency of the perturbation is much 
smaller than the temperature. (More accurate definition 
of the high temperature limit is given in Sec. V.) In this 
case the only contribution to the charge correlation func- 
tion, Q 2 , is given by the diagram shown in Fig. 3. The 
corresponding analytical expression is 
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Q2 = 4e 2 g dxdyK{x,y) d9 drF 2 (T 



4 



x V 



x + y 
2oj 



x + y t ^ x-y 



1 2lu 



x-y 

LO 



(26) 



2lo ' 2uj 
where 

lZ(x, y) — C\\ sin x sin y + C22 sin(a; + <j>) sin(y + 4>) 

+ C12 (sin(x + 4>) siny + sinxsin(y + tf>)) , (27) 

N\N T 



N, 



ch 



(28) 



is the dimensionless conductance through the dot from 
the left to right leads and 

Tt 



F(t) 



(29) 



sinh 2ttTt 

Expression ( ^6|) can be computed for different values of 
parameters. In Fig. 4 we present the result of computa- 
tion of Q 2 for two frequencies w — 0.l7 csc and lo = 7 CSC - 
Both those curves exhibit C 2 and \[C[ dependences at 
weak and strong pumping respectively. We also show 
the analytical curve given by Eq.(|37|) for the u> — ► 00 
limit. Below we discuss different limits of Eq.(|2q) 




-v- 


co =0.1 v 

'esc 


-A- 


'esc 




(t>»V 

'esc 



10 ' 



10 



C/y 



10 



10 



FIG. 4. The amplitude dependence of the pumped charge 
for different values of the frequency of the pump. For u) > 7 CSC 
the curves have the C 2 dependence at small values of C and 
the VC dependence at C > 7c SC . For small frequency (e.g. 
uj — 0.l7c SC ) there is an intermediate regime. 



A. Bilinear response 

First we consider the weak perturbation and perform 
an expansion of the diffusons up to terms linear in CV,-. 
As a result we obtain: 

/>oo 

Q2~ = 4ir 2 e 2 g d9e- 2N * e <1tF 2 {t) 
Jo J-e 



Cf{2uj0 - sin 2oj9) + C 2 sin 2uj8 



(30) 



where we have introduced the linear and circular pump- 
ing amplitudes 



C\ = Ch + 2C12 cos 
C c = 2sin^/c n C : 



a 



22. 



22 



°12- 



(31) 
(32) 



In the case of temperature T larger than the escape rate 
= Nc^dx/2%, (T ^> 7 CSC ) we find 



Q 24 5 T 7e 2 sc + 



^ 2 r 2 

Tesc 



C; 



(33) 



The second term of Eq.(|33J) survives the limit lo — > 0, 
thus reproducing the known result for adiabatic pump- 
ing [0J|]. On the other hand, this term vanishes at high 
frequency. The linear term is quadratic in frequency at 
small frequency and tends to a constant at large fre- 
quency. 

The linear pumping amplitude C\ in the case of two 
pumps has the form of Eq.(|3l]), which implies that the 
amplitude is just a vector sum of different pumps in the 
parameter space. On the other hand the circular ampli- 
tude is related to the area in the parameter space, covered 
by the pumps. 



B. Low frequencies 



Equation ( pq ) in the adiabatic limit lo — > is in agree- 
ment with the results of Ref. ||. Namely, this expres- 
sion gives the same asymptotic behavior for the limits of 
weak and strong pumping. To demonstrate this we con- 
sider the special case of the C matrix having the form 
Cn = C22 = C and C12 = 0. In this case we obtain 



-r I 2 5i 2C + 

wesc V 7csc \7csc T" ) ) 

) z = — eg . 

12 T 



\/7esc(7esc + 4C) 



(34) 



As temperature drops down to T ~ 7 CSC = A^ c h<5i/27r, the 
variance of the transmitted charge saturates to 



Q 2 = -e 2 g 



Si 



C 2 



* 7csc \/7 csc (7c S c + 4C) 3 



(35) 



The authors of Ref. || showed that for strong pump- 
ing the mean square transported charge is proportional 
to the length of the contour in the parameter space, and 
does not depend on the particular shape of the contour. 
Equations d34| ) and (35) support this statement, since 
for C 3> 7esc they reproduce a \fC dependence on the 
pumping amplitude. In the opposite case of weak pump- 
ing Eqs. (p3) and (|3^) give C 2 dependence in accordance 
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with jl],D. To understand the strong pumping depen- 
dence, we consider a loop in the parameter space. See 
Fig.(5). We notice (see Eq.@) and Refs. @J|, that adi- 
abatic pumping can be related to a contour integral in 
the parameter space. At sufficiently strong pumping the 
system at distant points of this space is uncorrelated and 
the total contribution to the pumped charge comes from 
the uncorrelated pieces of the loop, being proportional to 
their number. 





X2 










"(a) 






(b) 

FIG. 5. In figure (a) a loop is shown in the parameter 
plane. The grid divides the plane onto pieces, so that parts 
of the loop in the different pieces give uncorrelated contribu- 
tions to the transported charge. Figure (b) shows the loop in 
the phase space for strong pumping. In this case the loop can 
be divided onto pairs, and the pairs are not correlated. On 
the other hand parts of the loop of one pair are close to each 
other, so they are strongly correlated. 

In the limit of low frequency and zero C c (single pump) , 
the mean square fluctuation of the charge per cycle is 
quadratic in frequency. For weak pumping the ampli- 
tude of charge fluctuations is determined by Eq. (|3^) with 
4> = for arbitrary frequency u>. (For a single pump, C\ 
is the only parameter.) For strong pumping (but still 
lo 2 C\ < 7^ sc ) we find 



5767T 74 c T V7 esc 



3/ 2 



(36) 



We explain this dependence on the amplitude of the per- 
turbation in the next section. 



C. High frequencies 

In the limit of high frequencies, T 3> u) 3> 7eso the 
variance of the transmitted charge is given by 



\/7esc(7esc + 2C[) 



12 V7esc(7e S c + 2C70 



(37) 



In the limit of strong pumping this expression has the 
\[C\ asymptotic behavior. The curve Eq.(|3lj) is repre- 
sented in Fig. 4 by the solid line. 



IV. PHOTOVOLTAIC EFFECT AS PUMPING IN 
PHASE SPACE 

In this section we discuss the mechanism of charge 
transport by a single pump at finite frequencies (i.e. the 
irreversible contribution to the d.c. -current). We show 
the similarity with the mechanism of adiabatic charge 
pumping, discussed in In the adiabatic approxima- 
tion the system's motion is considered in a parameter 
space. For finite frequencies the parameter space has to 
be extended to phase space, which contains not only the 
perturbation parameters but also their time derivatives. 

According to Eq. ([)]), the transported charge for one 
period is determined by 



de 



2tt 



(38) 



e J dt J dt' J dr 

, „ / t + t' t\ „ + / t + t' t , , 
xTr^ U — + - St U — --|A 



We use the Wigner transform for the scattering matrix: 

S(t,t') = J S(e,(t + t')/2) e ^-*')|:. (39) 

We consider charge pumping at high temperature 
(T oj). In this case the integration over r is lim- 
ited by the inverse temperature 1/T. On the other hand, 
the scattering matrix S(t, t) in the Wigner representation 
varies slowly with respect to its time argument t. This 
allows us to expand the scattering matrices in Eq. ([38]) to 
linear order in r. Using the unitarity of the scattering 
matrix we finally obtain 



' = e 



dt 



de 



2tt cosh 2 e/2T 



(40) 



Tr < A 



as (e, t) 



dt 



<S+ (e,t)-S(e,t) 



as* (e, t) 



dt 



This equation was used by Brouwer in W . (See also jlT| ) . 
The scattering matrix in the Wigner representation is a 
function of the perturbation itself and its higher order 
derivatives with respect to time. (See Appendix C.) In 
the adiabatic approximation the derivatives are neglected 
as being small to higher orders in frequency. Beyond the 
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adiabatic approximation, we have to include the deriva- 
tives. 

We demonstrate that the analysis of Ref. can be 
applied to our case. We assume, that there is a single 
parameter tp(t). Then, following Brouwer, Ref. [GJ, we 
introduce a vector field 



Pi(e,t) 



ImTr { A ^lr 5t(M) }' (41) 



where Xi = d 1 'tp(t) / ' dt l and i is a non-negative integer. 

In these notations Eq.(^TJ) for the transported charge 
Q is given by 



de 



1 



2tt cosh 2 e/2T f-f 

' 1=0 



J2P l (e)dX l 



(42) 



The loop integral in the above equation can be rewritten 
as a surface integral using Stoke's theorem. We develop 
our analysis for the transported charge to the lowest or- 
der in frequency, so that the scattering matrix depends 
only on Xo = tp(t) and X\ = <f{t). According to Stoke's 
theorem for this two dimensional space, we obtain 

Q = e / * l2 X /tw1 / dtpdtpU.{e) , i 48 ! 



n(e) = Im Tr A 



2tt cosh 2 e/2T 

dS(e) dS^(e) 



dip dip 



In appendix C we present a formal derivation of OS / dtp 
from the equation of motion, Eq.(|l2|), in terms of the 
Green's functions of the dot to lowest order in oj/j csc . 

Now we interpret results found in the previous sec- 
tion using Eq.(|45j). For weak pumping we keep the time 
dependent perturbation to lowest order to calculate the 
derivatives with respect to ip and (p. Then we consider 
n(e) to be a constant and the integral over ip and tp 
gives the area of the contour in phase space. For a har- 
monic field tp(t) = coswi, the contour is an ellipse with 
large semiaxis y/C, small semiaxis wyC and area txujC . 
For the variance of the transmitted charge, we expect 
Q 2 cx w 2 C 2 /7 4 sc , which is in agreement with Eq.(|33|) for 

W < 7csc- 

In the limit of low frequency but strong pumping, we 
can apply Eq.(S) to understand Eq.(|36|). The power 
dependence [C^] is different from the adiabatic case 
[C 1 / 2 ]. The loop in phase plane is long along the ip axis 
but narrow in the tp direction because the frequency is 
small. [See Fig. 5(b).] The charge variation is determined 
by a sum of independent contributions from pieces of the 
contour along the tp axis. As can be seen from Fig. 5(b) 
the number of the independent pieces is iVj nc i = yj C /"f C sc- 
In the ip direction the system is correlated inside each 
piece of the contour since all points along the tp direction 
are separated by a distance, smaller than the correlation 
length. The characteristic area S c of each part is pro- 
portional to ujy/Cj csc . The variance of the transported 
charge can thus be estimated as 



Q 2 oc e 2 N ind S* cx e 2 



l 2 

/esc 



C 

Tcsc 



3/2 



(44) 



When the amplitude of the field C or the frequency 
w increases further, so that lo 2 C\ > ^ sc , this picture is 
no longer valid. The trajectory does not have parts close 
to each other and each part gives an independent con- 
tribution. The situation is similar to the case of strong 
adiabatic pumping, as shown in Fig. 5(a) and discussed 
in (|] . The variance of the transported charge is propor- 
tional to the total number of uncorrelated parts, so that 
Q 2 cx y/C, see Eq. @. 



V. LOW TEMPERATURE 

The previous discussion of d.c. -current generation is 
quite general. However, it does not take into account 
the heating of electrons by an external field which be- 
comes important at low temperature. In this regime, the 
electron distribution function in the dot changes and ac- 
quires a width larger than the electron temperature in 
the leads. 

The new width Th of the distribution function can be 
estimated from the following picture. An electron has 
random transitions between different energy levels. The 
time between consecutive transitions t tr is determined by 
the Fermi golden rule: 



2n\V nm \ S(e n - e m ± w) ~ \V nrn \ 2 /di 



C 



(45) 



The first equality sign follows from the Fermi golden 
rule, the second sign represents an estimate of the char- 



acteristic value of the matrix elements \V n 



and the 



density of states l/6t, the last equation is the defini- 
tion of C, cf. Eq.(pr|). Since an electron stays in the 
dot for a time r csc = 7"^ = 2iv/N c \ l 8i 1 it performs 
— lAtr7csc = C/7 CSC transitions. Each transition 
changes the energy of the electron by w. As in the ran- 
dom walk problem, the displacement of electrons in en- 
ergy space is cx u\fC/^ csc . 



This analysis gives a new temperature scale X^: 



C 



7e 



(46) 



This scale has a meaning only for strong fields, C 3> 7 CS c, 
so that the diffusion picture in energy space is valid. Oth- 
erwise, electrons experience few transitions with change 
of energy u>. Now we consider low temperatures, so that 
T ^ Th is not valid. We calculate the fluctuations of d.c. 
current for a system with a single pump. As we know 
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from Sec. Ill, at high frequency the number of pumps 
is not important and the result depends on their linear 
combination. 

Unlike the diagram, shown in Fig. 3, diagrams pre- 
sented in Fig. 6 have additional diffusons dressed on the 
distribution functions /(r). Collecting diagrams in Fig. 3 
and Fig. 6 we obtain the following expression for the vari- 
ance of the pumped charge: 



Q 2 = 4e 2 CW chff 



,11,11' I iW 



x V 
x V 



t + t' 
2 

t + t' 



t + t' 
2 

t + t' 



(1tF 2 (t) 
+ T,t-/j (47) 

- T, t - A 

r + co r+oo 

x / dgD(t,t-€,2r) / d£'V(t',t' -£',2t) 
Jo Jo 

x {2Csin 2 uo(t - f) sin 2 uo{t' - £') sin 2 lot 
+ N c h. sin cot sin tot } 

where g is the dimensionless conductance of the dot (see 
Eq.@) and 



F(t) 



T sin lot 
sinh 2ttTt 



(48) 



At high temperature T 3> w Eq.(|47]) reduces to 
Eq. (^6|) . Indeed, all three diagrams in Fig. 6 are smaller 
than the diagram in Fig. 3 at least by one factor lo/T. 

Now we discuss the limit of high frequency lo 3> 
max{7 csc ,C}. This inequality allows us to perform in- 
tegration over £ and £' in Eq.([l7|). We can replace 
sin 2 oj(t — by 1/2 and use the approximation 



/>oo 

Jo 



1 



Nrl 



2C sin lot 



(49) 



In this limit the product of the diffusons in the second 
and third lines of Eq.(p7|) does not depend on r. 

Energy of an electron in the dot changes due to the 
external field resulting in the redistribution of the elec- 
trons in the energy space. The new distribution func- 
tion becomes wider than that of electrons in the leads at 
temperature T. Consequently, the Fourier transform of 
the electron distribution function becomes narrower. The 



right hand side of Eq.(49) represent the effect of heating. 
This function appears in the integral over r in Eq.([47j) 
along with the function F(t), defined by Eq.(p9|). At suf- 
ficiently low temperature the convergence of the integral 
over r is determined by the 'heating factors', Eq.p9|) 
rather than by F(t). We note, that the shape of the 
new distribution function is not a Fermi function with a 
higher temperature. Instead, its Fourier transform has 





^ — +^^^ — 



F I ( ! . 6. Diagrams, which contribute to the d.c.-current at 
low temperature limit. (We do not show diagrams, which can 
be obtained from the above by omitting the upper or lower 
diffusons.) The last diagram contains the Hikami box, which 
is presented in Fig. 7. 

To be more specific, we consider the strong pumping 
limit C 3> 7cso when the electron distribution function 
is determined by the new scale Th, see Eq.(B^). We find 



16 LO 



(50) 



the form of the right hand side of Eq. (49) 



The same parameter dependence can be found from 
Eq. ( |37| ) replacing temperature T by the new energy scale 
T h . 

Equation (^7j) has two terms. One term contains fac- 
tor sin Lot sin Lot' . This term survives the high tempera- 
ture limit, see Sec. III. Nonetheless at low temperature 
the heating modifies the results of Sec. Ill, so that the 
temperature dependence saturates at the characteristic 
temperature scale Th. 



8 



The second term was completely neglected in the pre- 
vious sections. The origin of this term is similar to that 
of the thermoelectric effect in a conductor out of thermo- 
dynamic equilibrium. Although electrons are in equilib- 
rium in the leads, the heating changes their distribution 
function in the dot, producing a non-equilibrium distri- 
bution. Then non-equilibrium electrons escape from the 
dot. The direction of each escape is determined by the 
realization of the dot. An unbalance between electrons 
escaping through the left or right leads gives current. 
(sin 2 lot term in Eq. ( |47| ) reflects the electron-hole asym- 
metry, necessary for thermoelectric effects.) 




FIG. 7. The Hikami box, introduced in Fig. 6, can be 
obtain from these diagrams as well from their rotations. The 
grey rectangulars represent averages of the type GWG( R ) and 



described by a C 3 / 2 dependence of the variance of the 
d.c. -current on the field amplitude. 

We also considered a wide temperature range. At high 
temperature the variance of the current decreases as T _1 . 
We found that at low temperature heating becomes im- 
portant and it introduces a characteristic temperature 
Ih, see Eq. ([46|), below which the temperature depen- 
dence of the d.c. -current saturates. The result of heating 
on the d.c. -current is twofold. The first effect diminishes 
the d.c. -current by broadening of the electron distribu- 
tion function. The second effect produces the a ther- 
moelectric field, and is a non-equilibrium effect. This 
effect is related to the electron-hole asymmetry in the 
dot. Our results are different from the observed experi- 
mentally temperature dependence, see also Ref. ||. 

Finally, the photovoltaic effect is not symmetric with 
respect to inversion of magnetic field, similarly to adia- 
batic pumping case, see Ref. || . 

We are thankful to P. W. Brouwer for useful discus- 
sions. The work was supported by Cornell Center for 
Materials Research under NSF grant No. DMR-9632275 
(M.G.V. and V.A.) and Packard Fellowship (I.L.A.). 



VI. CONCLUSIONS 



APPENDIX A: 



We studied the d.c. -current through the quantum dot 
generated e.g. by time-dependent distortions of the dot 
shape. This d.c. -current is fluctuating from sample to 
sample and we found the second moment of its distri- 
bution. Unlike the previous works Jl],|3| on the adiabatic 
pumping we treated the system for a broad range of ex- 
ternal frequencies thus providing a bridge between adia- 
batic pumping and photovoltaic effects in microjunctions 
of Ref. §. 

The adiabatic approximation is not valid when the fre- 
quency of the perturbation to is comparable with the es- 
cape rate from the dot 7 0S c- Beyond the adiabatic ap- 
proximation the d.c. -current consists of reversible and ir- 
reversible parts. The reversible contribution was studied 
in the adiabatic regime. This contribution is determined 
by an integral over the contour in parameter space. On 
the other hand the irreversible contribution to the d.c- 
current is determined by an integral over this contour in 
phase space. 

A crossover from the bilinear (C 2 ) to \fC regime was 
found in ||. The crossover happens when the system 
makes a large loop during one period of the external field, 
so that it is uncorrelated at different points of the loop, 
see Fig. 5(a). We showed, that this crossover is universal 
and happens at arbitrary frequency. This result is consis- 
tent with the representation of the d.c. -current as an inte- 
gral along the contour in the phase space. An intermedi- 
ate regime exists for a single pump at low frequency and 
moderate amplitude of the external field. This regime is 



We define the wave function of electrons in channel a 
moving towards the dots by ip a (x,t) with x < 0, where 
|x| determines the distance from the dot boundary, see 
Fig. (8). Then tp a (x,t) for x > represents the outcom- 
ing electrons. The boundary x = is described by a 
superposition of the incoming and outcoming electron 
states and we denote it by tp a (0,t). The wave function 
of electrons in state i is denoted by ipi (t) . 



L 


s 




R 










incoming X<0 


— ► 











FIG. 8. Correspondence between the sign of x and the 
direction of motion of electrons with respect to the dot. 

We introduce the Keldysh Green's functions 

g a p(t,t',x,x')= (Ai) 

/ g™(t,if,x,x>) g ( a f(t,t',x,x')\ 
{ o g$(t,t,x,x') )' 



g\2\t,t',x') g%\t,t>,x>) 



(A2) 
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which are defined in terms of 

Gif (t, t', x, x') = -iQ(t - t')(ty a (x, t), i>\{x', t')}+), 
0$(t, t', x, x 1 ) = iG(t' - t)([^ a (x, i), ipl(x', t')}+), 

g ( a f(t, f , x, x') = -i([ip a (x, t), 4( x '> 01->> 

where [•, -]± denote commutator and anticommutator re- 
spectively. The similar expressions can be written down 
for gi a (t,t' ,x') Green's function, with ij) a {x,t) replaced 
by i>i{t). 

We assume that electrons do not interact in the reser- 
voirs and the Green's function of the incoming electrons 
(x; x' < 0) is given by the Keldysh structure: 



g a/3 (t,t',x,x') 

*/3 



(A3) 



G [ S{t-t',x-x') G%\t-t',x-x') 



V G^'(t-f,x-x') 

where 

G { *j{t,x)=iQ{t) 5 aP 8 (v F t - x) , 
G i^(M) =-*e(-*) 8 aP S(v F t-x), 
G { a f(e,x) = f a (e) (G${e,x) - G$(e,xj) , 



(A4) 
(A5) 
(A6) 



where /(e) is the distribution function of electrons in the 
channel a. In equilibrium with temperature T, 



f a (e) = tanh 



e - 8\i a 
2T 



(A7) 



where Sfj, a represent relative change in chemical potential 
for different leads. 

The equations of motion for the Green's functions de- 
fined by Eqs. (|Al[) and (A2) have the form: 



d 



d 



dt F dx 



g a p(t,t',x,x') (A8) 
S(x)W m g %l3 (t, t', x') + S(t - t')8(x - x% 

g ja (t,t',x') (A9) 



i| -*<,(*) 



= w^g Pa (t,t',o,x'). 

We notice that due to causality, G^g(t, t',0,x') = 
for x' < 0. This observation significantly simplifies fur- 
ther calculations. Indeed, we can represent the Keldysh 
component of the Green's function in the left hand side 



of Eq.(A9) in the form 



id/dt-H(t) 
x W^g af3 (h,t',0,x'). 



(Mi) (A10) 



The corresponding advance component is zero. Here 
l/(id/dt— H(t)) is the retarded component of the elec- 
tron Green's function in the dot. This definition is dif- 
ferent from that given in the main part of the paper, see 
Eq. (|l0|) . The latter will appear naturally in the end of 
this section. The additional term ~ W^W takes into 
account the escape from the dot through the leads. 
The next step is to represent Eq.(lAl) in the form 



'a/3 (t,t',x,x') - 
+ J dhdhG^it-t^x) 

xg^\t 2 ,t',o,x') 



g% ) {t,t',x,x')=G%\t-t',x-x') 

w- 1 



id/dt-H{t) 



(All) 

(*1,*2) 
7(5 



In the above equation we consider x = 0. Using 
G^g (t - t', 0) from Eq.© we find 



Gif(t,t',0,x')= Idh 



l-W- 



id/dt - H(t) 



x G ( s f(t 1 ,t , ,0,x'), x' <0. 



(Mi 



(A12) 



Substituting this expression to Eq.(All) and taking 
x = +\S\ — > 0, we obtain for x' < 0: 



S$?(tJ,+\6\y)= ldt 1 S cn (t ) t 1 )G$ ) {t 1 -1?,-x>), 

(A13) 



{K), 



where the scattering matrix S a p(t, t') is given by Eq.(|ll[). 

Equation (A13) is valid for x' < 0. We have to repeat 
the procedure described above to calculate the electron 
Green's function in the leads for x' > 0. Since the equa- 
tions which determine evolution of the Green's function 
from x' < to x' > are conjugated to those for x, we 
conclude, that 

glf(t,t\+\5l+\8\) (A14) 

f dhdhS^ (t, tl)G^P {h-h, 0)5jg (t 2 , t') 

The currents in the left and right leads are given by 
h{t) = ev F (A15) 

x^si^M+l^+l*!)-^*,*,-!*! 

Q = l 

I r (t) = ev F 

X 

a=JVi+l 



(A16) 



E (^(M +1*1, +1*1) - ^(M, -1*1, -|*|)J 



where 5 — > 0. This limit is just a reminder that 
Gaa(t,t, — \S\, — \S\) is taken for incoming electrons and 
is given by Eq.(A6). Consequently, 
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Qif(t,t-\s\,- 



/(+«>), 



(A17) 



/(*) = 



-HCxD 



Since the charge is conserved, I\{t) = — I r (i). We 
rewrite the current through the dot as 

7 W - = W ) = -m (A18) 



Subst ituting Eqs. ( |A14|) a nd ( |A17j ) into Eqs. ( |A15[ ) and 
( |AT^ ) and using Eq.(|A18|) we obtain Eq.(g). 



APPENDIX B: 

In this appendix we derive Eq.(|l4|) from the general 
Eq. (|l3|). The only assumption we are using here is that 
the distribution function of electrons is the same in all 
channels, i.e. f a (t) = f(t). Substituting the explicit 
form of the scattering matrix from Eq. (O) , we obtain 



0a 

(Bl) 

pa 



I{t) = 2mve!.Jdtif{t-t 1 )k ol p \ W^G {A) {h,t)W 
r dt 1 f(t 1 - t)A aP [w^G^(t,h)W 
2-niv [ [dt 1 dt 2 f(h - t 2 )A al3 

,t)W 

0a 



W^G (R] (t, h)WW^G {A) (t 2 ,t)W 



where the diagonal matrix A is defined in Eq. (|Tl|), and 
G(t,t') denotes the Green's function in the dot. 

Now we use the equation of motion of the retarded 
and advanced Green's functions, Eq. (|l^), to simplify 
the right hand side of Eq. (Bl). For this purpose, we 
pre-multiply the equation for the advanced component 
G^ A \t 2 ^t') by the retarded component & A \t,ti), then 
we post-multiply the transposed equation for the re- 
tarded component G^ A \t, h) by & A \t 2l t'). Subtracting 
from the second equation the first one, we obtain 



2mvG {R) {t, t x ) WW^G^ A \t 2 ,t') 
~d_ _d_~ 
dfi + dT 2 

+ & R) (t,h) [A(ti)-H(t 2 ) & A \t 2 ,t') 



(B2) 



& R \t,h)& A Ht 2 ,t>) 



& R) (t, h)5(t 2 - t') - S(t - h)G^ (t 2) If) 



Substituting this equation to Eq. 
Eq.ph. 



we obtain 



2eG(e,t)- Ho-mvWW\G{e,t) 



(CI) 



^ {2i) k k\ \ de k { ' de k J 



Here G(e, T) is the Green's function in the Wigner vari- 
ables e and t. We represented H(t) in the form H (t) = 
H + V<p(t). 

In the adiabatic limit only the k = is taken into 
account. This approximation is crucial in the case of a 
single pump. By appropriate choice of the beginning of 
the cycle the pump moves for the second half of the cy- 
cle along the same trajectory as for the first half, but 
in the opposite direction. As a consequence, the total 
transported charge Q, Eq. ([icf) , vanishes in the adiabatic 
approximation. To remove this symmetry, we can add 
another pump oscillating with the same frequency, but 
with different phase shift. Also, the higher order terms 
in Eq.(Cl) break this symmetry. 

We consider contribution to the lowest order in fre- 
quency to the transported charge Q. For this purpose, 
we neglect all terms with k > 2, keep the k = 1 term to 
the first order and include all orders in k = term. The 
solution is 



G (e,t) 



1 



e-H - V<p(t) + mvWW^ ' 



(C2) 



Gi(e, t) = i^- (Go(e, t)VG 2 (e, t) - G%(e, t)VG (e, t) 



To the lowest order in ip(t), the scattering matrix S 
has the form 



where 



and 



S 1 (e,t)=S (e,<p(t))+<p(t)A(<p(t)), 



S (e, <p(t)) = 1 - 2 7 n4F t G (e, t)W, 



j. . dS 



(C3) 



(C4) 



(C5) 



= - ( G (e, t)VG 2 (e, t) - G 2 (e, t)VG (e, t) 
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